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Extending the previous 2-gender dioecious biploid gene-mating evolution model, we attempt to 
answer “whether the Hardy-Weinberg global stability and the exact analytic dynamical solutions can 
be found in the generalized N-gender polyploid gene-mating system?” For a 2-gender gene-mating 
evolution model, a pair of male and female determines the trait of their offspring. Each of the pair 
contributes one inherited character, the allele, to combine into the genotype of their offspring. Hence, 
for an N-gender polypoid gene-mating model, each of N different genders contributes one allele to 
combine into the genotype of their offspring. We find the global stable solution as a manifold in 
the genotype frequency parameter space. We exactly sovle the analytic solution of N-gender mating 
governing equations and find no chaos. The 2-gender to N-gender gene-mating equation is analogues 
to the 2-body collision to the N-body collision Boltzmann equations with discretized distribution. 

PACS numbers: 87.10.-e, 87.23.-n, 83.80.Lz, 05.45.-a 


I. INTRODUCTION 

In a previous work[l], we had surveyed the time-evolution of genotype frequency in a 2-gender dioecious biploid 
population system. We derive our governing equations based on our four fundamental assumptions, 

(a) The mating of gene-holders is (approximately) closed in a population system. 

(b) The mating probability for certain gene-type to another gene-type, due to average-and-random mating mechanism, 
is proportional to the product of their population (quadratic form for 2-gender, cubic form for 3-gender, etc.). 

(c) The probability of a gene-type for newborn generation obeys the Mendelian inheritance. 

(d) The accumulation of human population obeys the exponential growth law, for both the birth and the death 
evolution. 

We had derived the governing equations in the parameter space of the genotype frequencies, and had solved the 
analytical solution of the governing equations within a single locus and arbitrary number alleles mating process. 
However, our previous work only focus on 2-gender mating, normally one male and one female. We can analytically 
prove the global stability and the Hardy-Weinberg equilibrium [2, 3]. For biological and mathematical interests, we 
ponder whether the N-gender case, where N is an arbitrary integer number, will bring interesting results. In other 
words, we like to consider a mating system of N different genders. The mating process of the system will be successful 
if and only if there is a group of N different genders participate in the mating process. For a certain genotype (or a 
biological trait), each of these different N genders contributes an allele, together combine into the inherited genotype 
of their offspring. 

The questions we ask and aim to address are: Whether 2-gender mating is special for N=2, so that it assures the 
global stability for the system under time evolution? Whether 2-gender mating is special so that there exists exact 
analytic solutions? Whether there is chaos for N-gender mating system at some integer N? Is it possible to solve the 
analytic solution for arbitrary N? Is there a specific integer number M, such that for gender number smaller than 
M, the system has the global stability; if the gender number is larger than M, then the chaos appears (from global 
stability to chaos)? If so, what is the integer M? 

Since I-gender mating is just trivial asexual reproduction, and 2-gender mating is a case with global stable solution, 
one might expect that 2-gender mating is a transition. In other words, for N> 2, there might have chaos or more 
complicated time-evolution phase diagrams. However, as we will show later, for the general N-gender gene-mating 
system, there is always a global stable manifold as long as the system satisfies the previous four assumptions. Through¬ 
out the paper, we focus on the standard equations without considering the mutation and natural selection. From 
the mathematic viewpoint, the highly symmetric governing equations result in a continuous global stable manifold. 
This high-dimensional curved manifold is a global stable attractor, under time evolution attracting all points in the 
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Euclidean fiber attached at every manifold point. The stable manifold and the fibers together form a fiber bundle 
completely fill in the whole genotype-frequncy parameter space. We solve the analytical solutions exactly for the 
general N-gender (n+l)-alleles gene-mating system. 

One remark is that the coefficients of the governing equations for N-gender mating can be derived explicitly from 
the well-known multinomial theorem: 
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( 1 ) 


Therefore, in the discussion below we will not explicitly write down the coefficients, except for the simple case of 
N = 1,2. Another remark is that our governing equation shall also be viewed as generalizing the biploid system to 
the N-polypoid system. For the N-gender-mating N-polypoid system, we believe that if the previous four assumptions 
are obeyed, the global stability of genotype frequencies is robust. 


II. 1-GENDER GENE-MATING MODEL (ASEXUAL) 

Here we start from the asexual mating system. There exists only one gender in a system, such that each gender 
could self-produce its offspring exactly the same as itself. Then we denote the population number of certain genotype 
as Gq,, where a is arbitrary interger number. Thus, parallel to the work [1], based on the assumptions done in Sec. I, 
we can derive the governing equations for population 

jG^= hG^-kdG^, (2) 

where kh is birth rate and kd is death rate. Apparently we can normalize the Gq, by the total population P, and 
we will relabel ^ —>■ Gq as the genotype frequency. Hereafter we will re-parameterize the equation in terms of the 
genotype frequency Gq, and we obtain the governing equations for genotype frequencies 

(3) 

Hence, any point in the genotype frequency parameter space is stable for 1-gender asexual system. The genotype 
frequency is within the range of [0,1]. 


III. 2-GENDER GENE-MATING MODEL 


This section briefly summarizes the analytic solution studied in our previous work[l]. We had studied the governing 
equations for genotype frequencies: 


^Gaa— kb{G^^ ^ Gaj + \ ^ ^ GaiGaj — Gaa), 
j—0 i—Oj—0 

n n n n 

-^Ga/3= kb{2G^^Gpp+ ^ GoiiGpp-\- ^ GqqG^j-I-^ ^ ^ GaiGpj—Gap)- 

i =0 j =0 z=0j=0 


( 4 ) 


Recall that the death rate kd has no net effect on the genotype frequencies, because it is a universal linear effect on 
each genotype population which does not change genotype frequency. Note that Gqq means the genotype frequency 
contributed from the same alleles a from both parents. Gap means the genotype frequency contributed from the 
alleles a and /3 from each of the parents respectively. Hereafter we will identify Gap and its Gpa as the same genotype 
frequency, therefore the governing equations have taken account of the combined effects of Gap and its Gpa- We had 
found an exact analytic parameterization of the global stable manifold: 

{ r' (f) 'i — r n sin^(^i) 1 cos^(6>„ + i) 
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Her 6i is within the range [0,7r/2]. More conveniently, we find that a new parametrization of coordinates, 


G aa — Gcta,eq (0,)+ E 

i—0,i^cic 

Ga.j^ = Goij3,eq 

providing an inverse transformation mapping between the two kinds of parametrizations: 

6»fc=tan“i{[X;(2G'--+X]Gij)]/(2G'j,_;^^_;^+ Y. Gk-i, j)} , 1 < A: < n, with n formulas 
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Sij can be solved from the following equations: 
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1 < fc < n, with n equations 


n 

tan (6>j,) = [Y (G„i+2smi)]/(G'^,fc_i+2sm,fc_i), 2<k<n, 0<m<k-2, with (”) equations 


i—k 


Given any given initial values of Gij parameters as the input data and the starting point of time-evolution, we are 
able to find the final equilibrium point solution Gij^eq- 
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By the above relation, we can simplify the inverse trnasformation, the function Saji{Gi 


e^+i=iim i{[ X) (2G„„-fEG'aO]/(2G +EG;xi)}i 0 </r < n - 1, with n formulas 
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Sa 0 — ^(Gajs^eqiGij) — Gafi) — |(|(2Gq,q, + Y G ai) {“^G p i3 + YG/Si) — Gafi), with formulas 

2=0 2=0 


(9) 


The importance of the new parametrization is that the original governing equations are simplified to a set of decou¬ 
pled governing equations in the new coordinates: 


r ^ -ktSa/ 3 , 0 < a < 13 <n 
1 ^=0, l</r<n 

Since the decoupled governing equations are the simple exponential-decay type differential equations, we obtain the 

exact solution of governing equations: 

n n n 

Gij{t)Gij= Gij^eq{diJ.=0ii)Gij+ SoipG {Gaa + G pp—2Gap) ■ (11) 

0<i<j<n 0<i<j<n O<a<0<n 

In the following sections, we will generalize the 2-gender gene-mating model to N-gender gene-mating model. 


IV. 3-GENDER GENE-MATING MODEL 


Hereafter we do not bother to write down the coefficients of the governing equations, which could be easily deter¬ 
mined by the multinomial theroem Eq.(l), the basic mating rules and the Mendelian inheritance. The governing 
equations of genotype frequencies for 3-gender mating have the form: 
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dt 
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dt 
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( 12 ) 


where ii < 12 ^ * 3 ! *4 ^ *5 ^ * 6 ! *7 ^ *8 ^ * 9 - We should remember that there will be a permutation factor for 
Gi-yigigGi^igigGi^i^ig, ^ geHotypos of throo genders are not the same. Again, we note that Gap-y means the genotype 
frequency contributed from the alleles a, /3 and 7 from each of the three genders respectively. Hereafter we will 
identify Gap-y and its index permutation Gpa-y, Gp^a, etc., as the same genotype frequency. In our notation, we will 
define the genotype frequency Gpo,-y such that a < /3 < 7 . 

Interestingly, regardless the non-linear governing equations are generalized to the more-complicated 3-gender cubic 
form (instead of the 2-gender quadratic form in Sec. Ill), we can still solve the governing equations and its global 
stability. We find an analytic parametrization of a continuous global stable manifold as follows: 
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For the later convenience to decouple the governing equations, we find that a ne-w parametrization of coordinates: 

71 n n 

Gcicict — Gcx.a.oL,eq (^/x) “1“ oli H“ Sgn(z O^) 'IL\ cx.i ^aij: 

^= 0 , 27 ^ 0 ; 2 = 0 , 27 ^ 0 ! 0 < 2 <j<n 

(14) 

GcaalS — Gaa^,eq (^/i) a/3 3ui oij3-i 

Gafd^ — Gcx^0^eq (^/x) H“ 3lii 

^ Gol^j — Got^-y^eq (^/x) 1 


where a < /5 < 7 . Given any initial input Gijk parameters, we are able to find the final equilibrium Gijk,e 
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• {2>G^^^ -1- 2 ^ Gppi -l- ^ Gjij) 

2=0 0<i<j 

i,j^a 

By the above relation, we can define the inverse transformation from the old to the new coordinates: 

72 72 72 72 72 

6 »^+i=tan-i {[ ^ (3G^„„+2 ^ G^ai + E G„y)]/(3G^pp + 2 E Gp^r -f E G^i^)}, 0 <^<n-l, n formulas 

a=/x+l 2=0 0<i<j i—0 0<i<j 

i^ct ijj^a i^P' 

•^1 2 ^^*^G(^,eq{,Gijk) Ga^/ 3 ,eg(G 2 j/j;) Gaa/3 Go^j^^')^ with ( 2 ) formulas 
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Tfiot^y — ^(Ga^-y^eg {Gijk) Gct^y)t with ( 3 ) formulas 


(16) 
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Our next key step is deriving the decoupled governing equations in the new coordinates: 


= -kbSi ap, 0 < a < P < n, 

= -kbUi a/3, 0 < a < P < n, 

= -kbmapj, 0<a<P<j<n, 

^ = 0 , 

Finally we obtain the exact analytic solution of governing equations for a 3-gender mating system: 


(17) 


^ ^ ^ ^t^')Gijk~k ^ ^ -Si a/ 3 C (Gaaa Gf^aji 0 a/ 3 / 3 “l“ 0 / 3 / 3 / 3 ) ( 18 ) 

0 <i<j<k<n 0 <i<j<k<n O<a< 0 <n 

n n 

-f ^ a/3^ (^acta 8(7aa/3“t“8t7a/3/3 ^/3/3/3)“t ^ ^ (^aaa“t“^/3/3/3 “t“ G^^^ 3Gfy^^P 

0<a</3<n 0<a</3<7<n 


Due to the exponential decay to the equilibrium manifold, we have proved the global stability of 3-gender mating 
system. 


V. N-GENDER GENE-MATING MODEL 


Now we will generalize the above analysis to the N-gender gene-mating model. Each genotype frequency parameter 
can be written as, Gai ■ ■ ■ aia 2 ■ ■ ■ 020:3 • • • 03 .Om • • • ctm > each is within the range [0,1]. We find the existence of 



the global stable fixed points as a continuous manifold and the parametrization of the global stable manifold 
as: 
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where we define oq = —1. All 9^ are within the range [ 0 , 7 r/ 2 ]. 

We extend the stable manifold through the following parametrized Euclidean fiber at each fixed point on the stable 
manifold. This procedure lead to our new parametrization of coordinates. For A-gender mating of (n-fl) alleles, 
we seperate the parametrization of fiber space to two cases, for N is an even integer and for N is an odd integer. 

• (1) For N G even, N = 2k: 

There are 3 types of vectors spaning the extended Euclidean bundle space. 

(a) The symmetry type for a o /3: 


Si afi Gq, ^ • • • a/3 • • • /3 ■ ■ ■ aP ■ ■ ■ P' '^here 1 < z < A: 


( 20 ) 


(b) Th anti-symmetry type for a gg P: 


Ml 0/3 = ... ^-(fc + l-f)G^... a/3 • • • p + ik+l-i)G^ ■■■ap-- ■ p-^a ■ ■ ■ aP ■ ■ ■ / 3 > where 1 < z < A:-l (21) 
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(c) The mixed of multiple-genes (mixed type number > 3) vector: 

‘^Oil * * * 0;i0;2 * * * Ct2. ~ T ^2Ga2---a2 T * * * T * * * G;iO;2 * * * 0^2.’ 

rim rii rim 

m 

where n = ni + 712 + • • • + rim = '^^rii (22) 


• (2) N € odd, iV = 2fc + 1: There are 3 types of vectors spaning the extended Euclidean bundle space, 

(a) The symmetry type for a o /3: 


Si ap Gq, ^ Gq, .. ,,, p Gq, ... aji ■ ■ ■ P • • • a/3 • • • /3 ’ '''^here 1 < z < A: 


(23) 


2fc + 2--i 


(b) The anti-symmetry for a O /3: 

Ui aji — Gq, ... (2fc-|-3—2 z)Gq, ... ... ^-|-(2fc-|-3—2 z)Gq, ... ... p—G^ ... where 1 < z < fc (24) 


(c) The vector of mixed multiple-genes (mixed type number> 3): 


rri(X\ • • • aia2 • • • a2. Olm * * * — ^iG^i-.-ai T ri2G T * • • -b rimGa^...a^ NGq^-^ . . . aia2 • • • 02. Olm ' ' ' Cln 


where IV = ni + 112 + ■ ■ ■ + rim = 


(25) 


For example, we denote y c(z) G q,q, «... p as (c(0), c (1), • • • , c(n)). 


N 
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i=0 

N-i y 

For N = 2, we have Si ajs = (1? —2,1). 

For N = 3, we have s'! ajs = (1, —Ij —1,1), ui a/s = (1, —3, 3, —1). 

For = 4, we have S 2 a /3 = (0,1,-2,1,0), si a /3 = (1,0,-2,0,1), zzi a /3 = (1,-2, 0, 2,-1). 

For iV = 5, we have S 2 a /3 = (0,1,-1,-1,1, 0), si a /3 = (1,0,-1,-1,0,1), M 2 a /3 = (0,1,-3, 3,-1,0), mi a /3 = 
(l,0,-5,5,0,-l). 

For iV = 6, we have S 3 a /3 = (0, 0,1, -2,1, 0, 0), S 2 a /3 = (0,1,0, -2, 0,1,0), si a /3 = (1,0,0, -2,0,0,1), M 2 a /3 = 
(0,1, -2,0,2, -1,0), Ml a /3 = (1,0, -3,0,3,0, -1). 

The parametrization of the genotype frequency space for the N-gender gene-mating of (n-|-l) alleles system is like 
the parametrization of the fiber bundle space including both the global stable base manifold and the fibers attached 
on the manifold. We obtain: 

n 
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Given any initial values of the genotype frequency G... parameters, we are able to find the final equilibrium G,,,^eq 

via the inverse transformation: 


G, 
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ni\n2\■ ■ ■ Um'- N 


(TV-2) ^ 


il,--- ,iN-l=^a 
, m n 

■ +iN-l)J2 Ga.-o.n 
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(27) 


j—1 k—0 


0 < 2 i<---< 2 fc<n 
^ 1»" ■ 


By the above relation, we can simplify the inverse trnasformation from the old to the new coordinates: 

22 TL 22 22 

0^+i=tan ^{[ ^ (A^Gq...q, + (A^ — 1) Gq-.-cij + (A^ — 2) Ga-.-aiji^ + • • • + GQiji 2 ...i„_j)]/ 
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( 28 ) 





Next we derive the decoupled governing equations in the new coordinates: 


^2^ = -hsi a/3, 0 < a < P < n 
= -hui a/3, 0 < a < /3 < n 

^ = -hm, where m = . Um-y Qm 

^^ = 0, l<A:<n." 


( 29 ) 


Parallel to Secs. Ill and IV’s analysis, the decoupled governing equations are the exponential-decay type differential 
equations, we obtain the exact analytic time-evolution solution of governing equations: 
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VI. CONCLUSIONS 

We have found the exact analytic solutions and proved the global stability of an N-gender gene-mating N-polypoid 
system. Since there is no chaotic behavior, it seems to suggest that the natural law does not go against N-gender 
N-polypoid system as far as the global stability concerns. Presumably the extraterrestrial aliens or intelligence may 
have N-gender N-polypoid system with N> 3 while the mating system still enjoys the global stability. It is likely that 
when the mutation and the natural selection process sets in, it will alter the governing equations to be less-symmetric, 
and thus possibly result in richer or more chaotic behaviors. It will be interesting to study the perturbation away 
from the symmetric governing equations, with the help of our exact analytic solutions. 

We can also compare the 2-gender to N-gender gene-mating equation to the Boltmann-like equation. The 2-gender 
to N-gender gene-mating equations are basically the 2-body collision to the N-body collision Boltzmann equations 
with discretized variable distribution instead of continuous one. It may be interesting to find a set of exact solvable 
time-dependent Boltmann equation by generalizing our discretized variables to continuous functional variables. Such 
a derived Boltmann-like equation should still provide the global stability and a continuous fixed-point global stable 
manifold in the infinite-dimensional continuous parameter space. All these will be left open for future directions. 
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